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METHOD FOR CALCULATING ALLOWABLE CREEP STRESS 
IN LINEARLY INCREASING STRESS ENVIRONMENT 
by Charles L. Whltmarsh, Jr. 

Lewis Research Center 

SUMMARY 

An analytical technique was developed to calculate the allowable design stress for a 
specified creep limit of a material under conditions of linear stress buildup and constant 
temperature. The analysis utilizes constant-stress creep data in the integration of a 
stress -dependent creep rate over the material history. Allowable design stress limits 
for 1 -percent creep of refractory metals in the temperature range 800° to 1600° K are 
increased by 25 to 57 percent when compared with calculations based on the constant end- 
of-life stress condition. 


INTRODUCTION 

An accurate assessment of the allowable stresses that cause high -temperature creep 
in materials is important for space power systems currently under study. A great body 
of experimental creep data has been obtained for materials of interest under conditions of 
constant stress and temperature. The Larson-Miller (ref. 1) correlation, for example, 
which plots allowable constant stress against a combined time and temperature parameter 
for a fixed value of creep, is widely used for interpolation among the data. However, con- 
ditions of variable stress and/or temperature may be encountered in practice, and few 
creep data for such conditions exist. The situation that prompted the present work was 
the calculation of cladding thickness for a cylindrical nuclear -reactor fuel element de- 
signed to contain fission gases being generated at a constant rate. Internal pressure and, 
therefore, stresses are not constant (although temperature can be reasonably assumed to 
be), but increase linearly from some initial low value. Assuming that the end-of-life 
stress is constant over the entire reactor lifetime is unduly conservative, and the calcu- 
lational results of cladding thicknesses and core volumes are greater than necessary. 



II 


This assumption may severely penalize high -temperature-reactor concepts intended for 
space power. 

If the creep rate of a material under stress is assumed to be independent of prior 
creep sustained at other conditions, then total creep can be accumulated. This idea is 
not new (ref. 2) and has been experimentally verified for several materials (refs. 3 to 5). 
Recently, this technique has been used for the design of radioisotope capsules where 
helium, a decay product, builds up internal pressure (ref. 6). 

The purpose of this analysis is to apply the concept of accumulation of total creep 
for variable -stress situations and to use as the basis for the analysis experimental data 
obtained under constant -stress conditions. 


ANALYSIS 
General Case 

The basic assumption used in the following analysis is that the creep rate of a mate- 
rial under stress is independent of prior creep sustained at other conditions and that 
total creep, therefore, can be accumulated, according to the relation 

f 1 de = f f e(cr, T)dr (1) 

Jq J o 

where e is the creep of a sample; e(o, T) is the creep rate of a sample and is in general 
a function of the stress a and the temperature T to which the sample is subjected; and 
and Tj represent the final fractional elongation and final time, respectively. 

If the creep rate is now assumed to be constant for the time period of an experiment 
at constant stress and temperature (second- stage creep only), then the creep rate e can 
be written as the ratio of some fractional elongation to the time tgfa, T) over which 
this elongation occurred. Thus 

€(<T, T)=-_!°- (2) 

t Q (a, T) 


To facilitate the translation of these data to other values of e and t, the time t f (a, T) 
required for the design creep to occur can be defined as 


t f (a, T) = 


e f 

k(a, T) 


( 3 ) 
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where the subscript f is used to differentiate between design conditions and the experi- 
mental test conditions in equation (2). Substitution of this relation into equation (1) and 
integrating the left side yield 


€ f = 


1 



( 4 ) 


To carry out the integration, t^o, T) is written as an explicit function of stress and tem- 
perature. A useful relation to accomplish this is the Larson-Miller (ref. 1) correlation. 
Under those conditions where the Larson-Miller correlation can be represented by a 
straight line on a semilog graph, it can be written as 

log a = mP + log a = mT(K + log tg) + log a (5) 


where 


m slope of equation 
P Larson-Miller parameter 
a ordinate a intercept 

T absolute temperature 
K constant dependent on material 
tg time of exposure 

Rearranging equation (5) into an explicit relation of time tp yields 


t 


0 



10 


-K 


and combining equations (2), (3), (4), and (6) yields 


(6) 
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( 7 ) 



Implicit in the derivation of equation (7) is the use of constant -stress - constant- 
temperature experimental data to describe the creep occurring under conditions of vari- 
able stress and temperature. However, this procedure is proper only if the basic 
assumption holds that creep rate is independent of prior creep. It should be noted also 
that, in the derivation of equation (7), no assumptions are required on the variation of 
stress and temperature over the operating lifetime of the sample. These variables often 
can be described as functions of the operating history of the sample; for example, 
a = cr(r) and T = T(r), which reduce the integrand to a function of time r. Depending on 
the circumstances, equation (7) can then be integrated, either numerically or in closed 
form, to determine the allowable end -of -life stress for a fixed core life or vice versa. 


Linearly Varying Stress and Constant Temperature 

A closed -form integration of equation (7) is obtainable for the relatively simple and 
interesting case for which the operating temperature remains constant over the sample 
history, while the stress varies linearly with time. This case is applied to the calcula- 
tion of creep in the cladding of a pin -type nuclear -reactor fuel element operating at con- 
stant power and constant temperature. Since fission -product gases are generated at a 
nearly constant rate, the pressure and, hence, stresses within the fuel element can be 
considered to build up linearly with time. Under this condition, the rate of increase of 
stress is constant and can be set equal to the stress at the end of the operating life di- 
vided by the time period over which this rate occurred. Thus, 

(8) 

dr Tj 


where the implicit assumption has been made that the stress is initially zero. 

Equation (8) can now be used to effect a change of variable of integration in equa- 
tion (7) from operating time t to stress a. Equation (7) then becomes 


a, = f° 10 K T f a 1/mT p cr-'/mT 


da 
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Integration yields 


|3 


Of = 



( 9 ) 


/ K \mT 

At this point, examination of equation (6) reveals that the quantity a (10 Tj J is equal 
to the stress value predicted by the Larson-Miller correlation for the condition £q = e^. 
This point is denoted as cr LM and a stress modifying factor ip is defined as 



( 10 ) 


Thus, equation (9) becomes 


a f “ ct LM^ 

and states that, when the stress is linearly increasing and the temperature is constant 
over the operating history, the allowable end-of-life stress to achieve a desired 
creep value is equal to the product of the stress cr^j^ (obtained from the Larson-Miller 
correlation for an equivalent value of creep) and a factor \p that depends on the material 
and the operating temperature. 

The \p factor was evaluated for several high-temperature refractory alloys of in- 
terest and is shown plotted as a function of temperature in figure 1. Data for these ma- 
terials were taken from references 7 and 8 and were assumed to have a constant Larson- 
Miller slope over the temperature range plotted. For the temperature range 800° to 
1600° K, the least sensitivity is exhibited by NAS-36 with ip values of 1. 25 to 1. 4, re- 
spectively, and the greatest sensitivity is exhibited by FS-85 with i p values of 1. 37 to 
1. 57, respectively. From figure 1 it can be seen that \p increases with operating tem- 
perature and is nearly linear over a broad range. 

To illustrate the use of the ip factor, imagine a hypothetical application requiring 
10 000 hours’ service at a constant temperature of 2200° F (1480° K) for which the stress 
on the sample is linearly increasing. If 1 -percent creep and the material T-lll are 
specified, for example, use of the Larson-Miller data given in reference 7 indicates a 
design stress limit a LM of about 5000 psi (0. 18 MN/m 2 ). From figure 1, ip for T-lll 
at 2200° F (1480° K) is about 1. 42. Hence, a stress cr f of about 7100 psi (0. 26 MN/m 2 ) 
at the end of life could be allowed. 
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Figure 1. - Calculational factor to increase allowable design stress for 1-percent creep in several refractory- metal 
alloys. 


Application to Reactor-Fuel-Clad Design 


In a study of a small, fast, gas-cooled, high-temperature reactor employing pin- 
type fuel elements, the cladding thickness was calculated for the constant stress 
as well as for the linearly varying stress. The allowable design stress and the pressure 
on the fuel cladding were related by means of the tangential stress equation for thick wall 
cylinders (ref. 9) evaluated at the inside surface of the tube wall 


ct = 



- 2r 2 P 
o o 





) 


( 12 ) 
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where r is the tube radius, P is the pressure, and the subscripts i and o refer to 
inside and outside conditions, respectively. Equation (11) was used to evaluate the 
allowable stress from which the cladding thickness was determined by equation (12). 

The cladding material was tungsten - 25 percent rhenium, and the operating conditions 
were 25 000 hours at 1800° K with 1 percent creep. The factor i p was 1. 6 under these 
conditions, with a 60-percent increase allowed in design stress over the constant- stress 
assumption. A saving in cladding thickness was effected, which allowed an increased 
fuel loading that resulted in a core volume reduction of 18 percent with the use of the 
technique reported herein as compared with that of the constant-stress assumption. 

The foregoing analysis utilized the following approximation to actual reactor condi- 
tions. Because the initial fission-gas pressure was zero and the gaseous coolant exerts 
a constant pressure on the outside of the element, the cladding is initially in compression. 
As the fission gases build up throughout the reactor life, this compression decreases and 
changes to tension when the fission-gas pressure exceeds the coolant pressure. At the 
end of life, the net pressure on the cladding equals the fission-gas pressure minus the 
coolant pressure. These pressure histories are presented in figure 2. The upper 
dashed line shows the fission-gas pressure alone, and the lower dashed line shows the 
net pressure (fission gas minus coolant) on the cladding. The solid line starting at zero 



Time 


Figure 2. - Pressure history of a nuclear reactor fuel 
element with respect to containment of fission gases. 
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and ending at the end-of-life net pressure P| was the assumed pressure history used for 
the stress calculations. This assumption is considered conservative, since the assumed 
tensile stress was always equal to or greater than the actual. 


DISCUSSION 

The analytical method presented herein appears mathematically suitable for any 
total creep limit within the constant- creep- rate assumption and is not restricted to the 
creep limit of the experimental test data. To assure such application, however, its use 
should be limited to low creep limits. Also, uncertainties exist in the use of the Larson- 
Miller creep rate correlation where longtime test data for materials of interest are 
usually nonexistent, and thus extended extrapolations are required. Fortunately, the 
stress-increase factor plotted in figure 1 is relatively insensitive to changes in both 
the value of the slope of the Larson-Miller correlation and the temperature. For exam- 
ple, a 100-percent increase in either m or T causes the factor to increase about 15 
percent. 

In theory, any stress-temperature- creep correlation and/or any stress history could 
be used if a closed analytical solution were not required. Thus, equation (4) could be 
solved numerically by use of a digital computer if sufficient test data and/or problem 
history data were available. 

Method application to nuclear-reactor design is restricted to those situations in 
which the fuel cladding thickness is determined by creep as opposed to those situations 
in which limitations are related to material elastic properties or strength. In general, 
this application will include high-power- density - high- temperature systems that are 
required to retain fission gases in the core. Of course, fission gas pressure is only one 
of the factors upon which cladding thickness depends, and in a detailed reactor design, 
additional items, such as thermal cycling, heat-transfer effects, and fabrication, must 
be considered. 
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SUMMARY OF RESULTS 

An analytical method was developed to calculate allowable end-of-life stress for a 
design creep limit of a material under conditions of linearly increasing stress and con- 
stant temperature. For a given material and creep limit, a temperature dependent fac- 
tor was calculated which represents the ratio of the allowable end-of-life stress to the 
constant stress obtained from Larson-Miller data. 

Calculations were performed to show that the allowable end-of-life stress for 
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refractory metals, for the conditions of 1-percent creep, and linearly increasing stress, 
and a temperature range of 800° to 1600° K, was increased by 25 to 57 percent when 
compared with constant-stress calculations. 

Also, the core volume of a small, fast- spectrum, high-temperature, gas- cooled 
reactor was reduced by 18 percent when the method presented herein was used for the 
determination of the required fuel cladding thickness. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 7, 1967, 

120-27-06-05-22. 
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